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ABSTRACT 

Infrared (IR) luminosity is fundamental to understanding the cosmic star formation history 
and AGN evolution, since their most intense stages are often obscured by dust. However, lo- 
cal IR luminosity function estimates today are still based on the IRAS survey in the 1980s, 
with wavelength coverage only up to 100/im. The AKARI IR space telescope performed all 
sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160/im) with 3-10 times better sensitiv- 
ity, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined 
with a better spatial resolution, AKARI can much more precisely measure the total infrared 
luminosity (Ltir) of individual galaxies, and thus, the total infrared luminosity density in the 
local Universe. 

By fitting modern IR SED models, we have re-measured Ltir of the IRAS Revised 
Bright Galaxy Sample, which is a complete sample of local galaxies with ^eo/im > 5. 24 J?/. 

We present mid-IR monochromatic luminosity (i^Liy) to Ltir correlations for Spitzer 
8/im, AKARI 9/im, IRAS 12/im, WISE 12/im, ISO 15/im, AKARI 18/im, WISE 22/im, 
and Spitzer 24/im filters. These measures of Lmir are well correlated with Ltir, with scatter 
ranging 13-44%. The best-fit Lmir-^o-Ltir conversions provide us with estimates of Ltir 
using only a single MIR band, in which several deep all sky surveys are becoming available 
such as AKARI MIR and WISE. 

Although we found some overestimates of Ltir by IRAS due to contaminating cir- 
rus/sources, the resulting AKARI IR luminosity function (LF) agrees well with that from 
the IRAS. We integrate the LF weighted by Ltir to obtain a cosmic IR luminosity density of 
nTiR= (8.5^2:3) X 10^ LoMpc-^ of which 7±1% is produced by LIRGs (Ltir > lO^^L©), 
and only 0.4±0.1% is from ULIRGs (Ltir > lO^^L©) in the local Universe, in a stark con- 
trast to high-redshift results. 

We separate the contributions from AGN and star-forming galaxies (SFG). SFG IR 
LF shows a steep decline at the bright-end. Combined with high-redshift results from the 
AKARI NEP deep survey, these data show a strong evolution of ^^ff^ (x(l+z)^-^^^-^, and 
^T?R oc(Uz)^-^^^-^. For n^f^, the ULIRG contribution exceeds that from LIRG already 
by z~ 1 . A rapid evolution in both ^t?r ^rm suggests the correlation between star 
formation and black hole accretion rate continues up to higher redshifts. We compare the 
evolution of ^t?r ^^^^ X-ray luminosity density. The ^t?r I ^x^^ay ^^^^^ shows a 
possible increase at 2: > 1, suggesting an increase of obscured AGN at 2: > 1. 
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1 INTRODUCTION 



To understand the cosmic history of star formation and black 
hole growth, we must understand infrared (IR) emission; the 
more intense star formation, the more deeply it is embedded in 
the dust, hence, such star formation is not visible in UV but in 
the infrared. Similarly, AGN evolutionary scenarios predict that 
they are heavily obscured at th eir youngest, Compton-thick stage 
(iTreister. Urrv. & Viranil [2009h . The Spitzer and AKARI satel- 
lites revealed large amounts of infrared emission in the high- 
redshift Univer se, showing stro ng evolution in t he infrared lumi- 
nosity density (Le Floe' h etaP 2005; Perez-Gon zalez et al ] |2005l : 
iBabbedge et al.l2006l:ICaputi et al. 2007; MagnelU et al.'2009'). For 
example atz=l jGoto et al.l (12010) estimated 90% of star formation 
activity is hidden by dust. However, the key baseline of these evolu- 
tion studies is still at z=0, using a local infrared LF from the IRAS 
survey in 1980s. 

For more than 25 years, bolometric infrared luminosities 
(i^T/i?,8-iooo/^m) of local galax ics ha ve been estimated using a 
simple polynomial presented by IPeraultl (119871) , obtained assum- 
ing a simple blackbody and dust emissivity. Furthermore, the red- 
dest filter of IRAS is lOO/xm, which does not span the peak of 
the dust emission for most galaxies, leaving a great deal of uncer- 
tainty. A number of studies f ound cold dust that can not be detected 
with the IRAS. For example, IPunne &Ealesl(l200ll) detected such 
cold dust with T y 2QK using SCUBA 450,850 /xm observation. 
ISvmeonidis et al detected cold galaxies with SED peaks at 

longer wavelengths using Spitzer/MIPS. These results cast further 
doubt on Ltir estimation made with only < lOO/xm photometry. 
More precise estimate of local Ltir and thus the local IR lumi- 
nosity function (LF) have been long awaited, to be better compared 
with high redshift work. 

' 1 

^AKARI, the Japanese infrared satellite (Murakami et al. 

l2007h . provides the first chance to rectify the situation since IRAS; 
AKARI performed an all- sky survey in two mid-infrared (cen- 
tered on 9 and 18 /xm) and four far-infrared bands (65,90, 140, 
and 160/im). Its 140 and 160/im are especially important to cover 
across the peak of the dust emission, allowing us to accurately mea- 
sure the Rayleigh- Jeans tail of the IR emission. The time is ripe in 
terms of modeling as well; se veral sophisticated infrared SED mod- 
els have become available (Charv & Elbaz' "2001"; 'Dale & Helou" 
20021 : iLagache. Dole. & Pu get 2003; Siebenmorgen & Krugel 
2007), from which we can not only obtain accurate Ltir esti- 
mates, but also constrain the astrophysics. 

In this work, we aim to re-measure local Ltir, and 
thereby the IR LF of the Revised Bright Galaxy Sample (RBGS, 
ISanders et al. 2003), which is a complete sample of local IR galax- 
ies. This work provides us an important local benchmark to base 
future evolution studies at high redshift both by the current AKARI 
and Spitzer satellites, and by next-generation IR satellites such as 
Herschel, WISE, JWST, and SPICA. We adopt a cosmology with 
{h, Qrn, ^a) = (0.75, 0.3, 0.7) for a comparison purpose. 
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Figure 1. Angular separation between AKARI PIS and IRAS sources. We 
matched sources within 60 arcsec of separation. 



over its precursor (IRAS), the presence of 140 and 160/im bands is 
crucial to measure the peak of the dust emission in the FIR wave- 
length, and thus the Ltir of galaxies, especially ones with lower 
dust temperature. 

We have cross-correlated the AKARI bright s ource catalog 
with the Revised Bright Galaxy Sample (RBGS; ISanders et alJ 
l2003h . using a matching radius of 60 arcsec, as shown in Fig[T] 
The RBGS is a complete, flux-limited survey of all extragalactic 
objects with 60/im flux density greater than 5.24Jy, covering the 
entire sky at |6| > 5 deg. All 629 objects in this catalog have mea- 
sured spectroscopic redshifts either from optical or millimeter/radio 
(CO or HI data). The maximum redshift of our sample is 0.087. 
Because of the completeness and availability of the spectroscopic 
redshifts, this catalog is suitable to construct local infrared lumi- 
nosity functions (LPs). In addition, each galaxy has a measured IR 
luminosity from the IRAS survey, allowing us to compare with the 
AKARI-based one. Due to the difference in sky coverage (IRAS 
covers 96% of the sky, and AKARI covers 94%), 24 galaxies out 
of 629 RBGS galaxies did not have an AKARI counterpart. We re- 
moved these galaxies from our analysis, but applied corresponding 
completeness corrections to our statistical analysis. 

One important caveat is that the RBGS is a flux-limited sample 
from IRAS (S'eo/.m > 5.24 J?/). Therefore, if the IRAS S'eo/.m flux 
was an overestimate as we discuss later, our work does not address 
incompleteness in the sample selection. However, at this very low- 
redshift (average redshift of 0.0082), the RBGS is still the largest 
complete sample of IR galaxies with spectroscopic redshift. 



2 THE AKARI ALL SKY SURVEY AND RBGS 

In this work, we use the /3\ version of the AKARI/IRC bright 
source catalog and the /32 version of the AKARI/FIS bright 
source catalog. The 5 a sensitivities in the AKARI IR filters 
(S9W, LISW, N60, WS, WL and A^IGQ) are 0.05, 0.09, 2.4 , 
0.55, 1.4, and 6.3 Jy dlshihara et al.]l201Ql : lYamamura et al.ll2009l) . 
In addition to the much improved sensitivity and spatial resolution 



3 INFRARED LUMINOSITY 
3.1 Estimating total IR luminosity 

For these galaxies, we estimated new total IR luminosities (Ltir) 
based on the AKARI photometry using the LePhare codeQ to 

^ http://www.cfht.hawaii.edu/~ arnouts/lephare.html 
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Figure 2. An exampl e of the SEP fit . We fi t the AKARI 6-band photometry 
to the SED model of .Charv & Elbad <200lh to estimate Ltir. 



fit the infrared part (>7/im) of the SED and estimate TIR lumi- 
nosit y. We fit our AK ARI FIR photometry with the SED templates 
fromlchary & Elbad (CHEL hereafter 2001). Although the shapes 
of these SEDs are luminosity-dependent, the large baseline from 
AKARI observations {S^W, LISW, A^60, WS, WL and A^IGO) 
allows us to adopt a free scaling to get the best SED fit, which 
is then rescaled to derive Ltir. Since the AKARI catalogs are the 
/3- version, we adopted a minimum error of 25%. FigO shows an 
example of the SED fit. At the median redshift of 0.0082 of the 
RBGS, peculiar velocity is not negligible. Therefore, instead of the 
measured redshift, we used the distance, which is corrected for he- 
liocentric redshift using the cosmic attractor model ( Mou ld et al.l 
I2OOO) for the SED fitting and to compute Ltir. In this work, Ltir 
is measured in the wavelength range of S-lOOO/xm. 

We chose not to use IRAS 12,25,60 and 100 /xm fluxes in the 
SED fitting. Due to the detector pixel size, AKARI has significantly 
better spatial resolution (30-40") than IRAS (1-2'). In a few % of 
cases, AKARI measured significantly lower flux than the IRAS, be- 
cause AKARI can subtract the background cirr us better, and/or can 
separate nearby confusing sources better. (See IJeong et all bOOVh 
for more details.) 

In Fig|3] we compare Ltir measured by AKARI with those 
by IRAS. The IRAS-based Ltir is measured using the following 
equation. 

Ltir{Lq) = 4.93 X 10"^^[13.48L^(12/im) + 5.16L^(25/im) 
+2.58L^(60/im) + L ^ {100 fim)]{ergs~^Hz~ I}) 

This equation is obtained by fitting a single temperature dust emis- 
sivitly model (e oc z^~^) to the flux in all four IRAS bands and 
should b e accurate to zb5% for dust temperatures in the range of 
25-65K (|Peraultl[T987b . The Ltir measured by AKARI agrees 
well with those by IRAS. However, the scatter increases toward 
lower luminosity in Ltir measured by AKARI than by IRAS, es- 
pecially at logiLTm) < 10. This is primarily because AKARI's 
higher spatial resolution can separate background cirrus or confus- 
ing sources better than IRAS, resulting in lower Ltir. 
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Figure 3. Ltir measured by the AKARI is compared with those mea- 
sured by the IRAS (Eq^ for the RBGS. 



Table 1. Co mparison of Ltir. e stimates with different SED models 
jCharv & Elb az 2001; D ale & Helo u 2002; Lagache, Dole, & Puget 20031) . 
FiglUpresents corresponding plots. 



Models 


cr(%) 


Offset (%) 


CHEL vs IRAS 


44 


-23 


Dale vs Lagache 


24 


24 


Dale vs CHEL 


10 


11 


Lagache vs CHEL 


22 


-13 



3.2 Model-to-model variation 

In this section, we check the SED model dependence of L tir 
measurements. The IR SED models we used are JCharv & ElQ 
(l200lh : lDale & Heloul (l2002h : lLagache. Dole. & Puget (2003). We 
used these models to compute Ltir exactly as described in Sec- 
tion l3.1l In FiglH we plot the measured Ltir against the model re- 
sults, and find a tight correlation. As shown in Table[T] the standard 
deviations between different Ltirs are 10-24%. The agreements 
are quite good considering the large flux errors 25%) associ- 
ated with infrared photometry. Median offsets between them are 
— 13 ~ +24%, which could also be caused by calibration prob- 
lems. In Table [T] the CHEL model shows smallest offsets from 
other models. For k-correction and for Ltir estimation (Section 
13. IK we adopt CHEL models, which give luminosities in between 
the Dale models (11% higher) and Lagache models (13% lower) . 

Overall, there are good agreements between models, suggest- 
ing that once you have FIR photometry (WL and N160 in this 
work), the model-dependent uncertainty in Ltir is small. This 
justifies our use of a single SED library for the final Ltir mea- 
surement. 

3.3 MIR-Ltir correlation 

A good correlation between Mid-IR monochromatic luminosity 
and Ltir is known to exist. This is especially true for SF galaxies 
because the rest-frame MIR luminosity is dominated by prominent 
PAH features such as at 6.2, 7.7 and 8.6 /xm, and the Ltir is by 
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Figure 4. Comparing Ltir measured using three different SED models dCharv & Elba3l 200 ll: Bale & Heloull2002l:lLagache. Dole. & Pugetlbooi) . Scatter 
and offsets between the models are summarized in Table[T] 



dust emission heated by the SF activity; i.e., both of these are good 
indicators of SF activity. Observationally MIR detectors are more 
sensitive than in the FIR. Therefore, if we obtain accurate conver- 
sions from MIR to Ltir, they will bring Ltir measurement to 
much more galaxies. 

AKARI covers both the MIR and FIR range with its 6 band- 
passes, measuring both MIR and total luminosity well. In addition, 
because it is an all- sky survey, it provides the 6-band photometry 
for ^600 RBGS galaxies. Here we have one of the best opportuni- 
ties to investigate the Lmir-Ltir correlation. 

We compute the MIR monochromatic luminosity at 8, 9, 12, 
15, 18, 22, and 24 /xm using the AKARI photometry at 9 and 
18/im, and the best-fit CHEL model from Section [TTI For the def- 
initions of monochromatic luminosity, we tried to choose luminos- 
ity in passbands of existing IR satellites, for precise and easily- 
reproducible definition, at the same time avoiding narrow emis- 
sion/absorption features adding too much noise to the resulting re- 
lation. In addition, this choice allows users to apply the relations 
easily to observed fluxes, at least for low-redshift galaxies whose 
k-correction is small. We checked this choice of monochromatic 
filters did not have much effect compared to the exact luminos- 
ity at a certain single wavelength. We later show that MIR-Ltir 
correlations for WISE12/j.m and IRAS12/j.m are almost iden- 
tical. Our choices of filters are Spitzer 8/im, AKARI 9/im, IRAS 



12/im, WISE 12/im, ISO 15/im, AKARI 18/im, WISE 22/im, 
and Spitzer 24/i m. 

For Lspitzer8fim, LAKARIQ^im, LiRAS12^im and 

LwiSEi2fim, we Started from the observed 9/j.m flux, then 
used the best-fit SED (to the 6 bands) from the CHEL model 
(§ 13.11 1 to color-correct the observed AKARI 9/im flux to 
Lmir in each filter. To be specific, we computed LspitzerSfim, 
LAKARi9fim, LiRAsi2fj,m and LwiSEi2^im through the filter 
response function of the Spitzer/IRAC S/xm, the AKARI 9/im, 
the IRAS 12/im, and WISE 12/im filters. Here, the best-fit model 
is only used to color-correct the observed 9/im flux, and thus, the 
obtained Lmir is not an integration of flux of the best-fit model. 

Similarly, for LIsol5^lm, LAKARiiSfxm, LwiSE22^xm, and 
Lspitzer24fAm, wc Started from the observed AKARI 18/im flux, 
color-corrected the flux using the color-correction from the best- 
fit SED model in § 13.11 then converted them to luminosity in each 
passbands. This color-correction is very small (^1% at most) be- 
cause AKARI's 9 and 18 /xm are very close to other passbands. 
Thus, the results in Fig|5]show almost observed mid-IR luminosity 
against observed Ltir- 

In FigJS] we show monochromatic luminosity in 8,9,12,15,18, 
22 and 24 /xm against Ltir measured using the CHEL model. In 
Fig|5] there exist good correlations in each panel. In fact, the cor- 
relations look like a simple scaling relation in all panels, i.e., the 
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stronger the MIR-emission, the larger the Ltir. The results sug- 
gest indeed monochromatic luminosity in mid-IR range represents 
the Ltir well. 

We fit a linear equation to the relation obtaining the following 
results. Note that the unit is in solar luminosity. 

The best-fit relation between Ls^m and Ltir is 



(20±5) xi/L°jt™-°'(±44%). 



Ltir 

The best-fit relation between Li2fj,m and Ltir is 
Ltir = (17 ±4) x i^L";™' (±25%). 

The best-fit relation between Li^^rn and Ltir is 



Ltir = (24 di 8) x iyL\ 



i/,15/j,m 



'(±13%). 



(2) 



(3) 



(4) 



(5) 



The best-fit relation between L2Anm and Ltir is 
Ltir = (51 ± 19) x i.L°J^±°-"(±24%). 

The logarithmic slopes of the relations are very similar to each 
other, perhaps reflecting the fact that most of these IR filters are 
wide and have some overlaps in wavelength coverage. This also 
means in this wavelength range (8-24/im), and within broad filters, 
the SED shape does not depend much on Ltir. However this re- 
sult is at low redshift, and whether this applies at higher redshift 
needs to be examined separately. We note that Liy,i5fim shows a 
less dispersed r elation with Ltir. 

For L ^n^ .lCaputi et all bOOVh (ealgl). iBavouzet et all (l2008h 
(eqEl), and lBoquien et al.l ( 201Ql) (eqlSI) present a correlation. 



Ltir = 1.91 x {uLjy) 



restSfim 



(±55%) 



Ltir = 377.9 x (z.L.)°f/,8/.m(±37%) 



(6) 



(7) 



(8) 



Ltir = 216.9 x (z.L.)?^2^8m^(±29%) 

Note that for iBoquien et al.l (l20ld) . we used their conversion in 
Table 1 from SINGS data. The unit is also converted from W to 
Lq. We overplot these relations in Fig|5] There is good agreement 
between these conversions and our best-fit relation, especially at 
9 < IoqLtir < 11, w here both sa n iples have enough num- 
ber of galaxies. N ote that ICaputi et al.l (l2007h 's sample spans at 
< ^ < 0.6 and lBavouzet et al.l fcOO sS's is mostly at z ^ 0.4. 
Thus, these samples are significantly at higher reds hifts than ours. 
SEP e volution could cause a small difference. See lBoquien et al.l 

for metallicity dependence of the conversion^ 

For Li2/2m, there exists a conversion from iTakeuchi et al.l 

(l2005h . 



log Ltir = 1.02 + 0.972 log Li2^,n 



(9) 



We overplot the relation in Fig|5] The slope is in good agreement. 
There is a ^50% offset in zero points. Possible reasons of this off- 
sets include different definition in Li2^im and different estimates 
for Ltir. Our Li2um rel ation is virtually identical to the one given 
bv lSpinodio et al ] ([T995l) which is, in our units: 



logLTIR = 1.51 +0.942 log Li2;.m, 



(10) 



For L2A um, we compare the conversions with the fo llowing in 
the literature jBoauien et al. 11201 d : lBavouzet et allboOSh . 

Ltir = 57.9 x (z^L.)^4Mm(±54%) (11) 

Ltir = 6856 x (z^L,)§4^^'^(=b54%) (12) 
Our conversion is in very good agreement with iBoquien et al.l 



(l2010h . I Bavouzet et al ] (l2008h 's relation has a slightly shallower 
slope, but in good agreement at 9 < IoqLtir < 11, where most 
of their galaxies lie^ 

Smith et al. (l2007h observed IT/xm PAH complex features in 
Spitzer IRS spectra of nearby star-forming galaxies. This feature 
was not detected previously with ISO, and therefore not i ncluded 
in the SED model we used. Although [Smith et all fcOOVh 's SED 
templates only contained two data points in the far-IR, and thus 
were not suitable for our purpose, we try to estimate how the 
17/im PAH complex features affect the Li^^m-LTiR conversion 
as follows. First, we remove a continuum from both of the spectra, 
jCharv & E lbaz 2001; Smith et al. 2007). Then we replace the 6- 
20jtxni region of the Charv & Elbaz (2001)'s spectra by Smith et all 
(l2007h . Here we scaled the spectrum so that the amplitude of the 
llfim PAH complex matches. Thus, t his ne w templates keep the 
overall SED shape of Chary & Elbaz ("2001"), but have the 17/im 
PAH complex features from Smith et al. (2007). We rerun the SED 
fit to estimate the Li^^m-LTiR conversion to find that Ltir es- 
timate (as a function of Lis^rn) will be ^3% smaller if the 17/im 
PAH complex features are included in the SED. Although the effect 
is smaller than the errors of the conversion, this indicates that future 
SED models need to include the 17/im PAH complex features. 

With current technology, detectors are much more sensitive in 
the mid-IR range than in the far-IR range. Therefore, once these 
correlations are assumed, they will provide a useful conversion to 
compute luminosities of many more galaxies, either at higher red- 
shift or at a fainter flux level, from a single MIR flux. For example, 
the above relations can provide a total IR flux measurement for all 
sources from the AKARI MIR all sky survey, which has 3 times 
more sources than the FIR all sky survey. WISE will perform even 
deeper all sky survey in 3.4, 4.6, 12 and 22 /xm in the near future 
( Wright 2008). For future use, we also provide Li2^im-LTiR, and 
L22^m-LTiR conversions in WISE filters as follows, and also in 
FigEl 

The best-fit relation between LwiSEi2iJ.m and Ltir is 
Ltir = (18 ± 5) x i/L°;^^±°^l2^^(±23%). (13) 

The best-fit relation between LwiSE22fj,m and Ltir is 
Ltir = (53 ± 20) x lyL'^H"^^'' (±25%). (14) 



4 INFRARED LUMINOSITY FUNCTIONS 
4.1 The 1/Knax method 

With accurately measured Ltir, we are ready to construct IR 
LFs. Since the RBGS is a flux-limited survey, we need to correct 
for a vo lume effect to compute LFs. We used the l/Vmax method 
(ISchmidtll9 68) for this. An advantage of the l/Vmax method is that 
it allows us to compute a LF directly from data, with no parameter 
dependence or a model assumption. A drawback is that it assumes 
a homogeneous galaxy distribution and thus is vulnerable to local 
over-/under-densities (Takeuchi, Yoshikawa, & Ishii 2000). 

A comoving volume associated to any source of a given lu- 
minosity is defined as Kiax = V^^max — ^^min' where Zmin is the 
lower limit of the redshift and 2;max is the maximum redshift at 
which the object could be seen given the flux limit of the survey. 
In this work, we set Zmin =0.0004 since at a very small redshift, an 
error in redshift measurement is dominated by a peculiar motion, 
and thus, Ltir also has a large error. This only removes 5 galaxies 
from the sample. 

For the RBGS, the detection limit is IRAS Seof^m =5.24Jy. 
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Figure 5. Correlations between Ltir and monochromatic luminosities (8,9,12,15,18, 22and 24 /im). The lines show the best-fit polynomial, based on w hich 
we present conversions. When available, conversions from the literature iTakeuchi et al...20Q5. : Caputi et al.. 2007 : Bavouzet et al. 2008 : Boquien et al...2Ql(t) 
are overplotted for a comparison. 
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We used the SED templates dCharv & Elbazll200lh for k-correction 
to obtain the maximum observable redshift from the flux limit. 
For each luminosity bin then, the LF is derived as 



E 



(15) 



t ' 



,where Vmax is a comoving volume over which the i\h galaxy could 
be observed, and AL is the size of the luminosity bin (0.3 dex). The 
RBGS is a complete in 60/im at IRAS Sm^im >5.24Jy. Complete- 
ness correction in terms of sky coverage of both satellites is taken 
into account. 



4.2 Monte Carlo simulation 

Uncertainties in the LF values stem from various factors such as the 
finite numbers of sources in each luminosity bin, the k-correction 
uncertainties, and the flux errors. To compute these errors we per- 
formed Monte Carlo simulations by creating 1000 simulated cata- 
logs, where each catalog contains the same number of sources, but 
we assign each source new fluxes following a Gaussian distribution 
centered at fluxes with a width of a measured error. Then we mea- 
sured errors of each bin of the LF based from the variation in the 
1000 simulations. These estimated errors are added in quadrature 
to the Poisson errors in each LF bin. 



4.3 IR luminosity function 

In Fig|6l we show LF of the RBGS but using AKARI photome- 
try. FiglT] shows number of galaxies used to compute the LF. The 
original LF measured using the IRAS photometry is also overplot- 
ted. There is a very good agreement between the LFs measured by 
AKARI and IRAS over the luminosity range of %<\ogLTiR <12. 
Although the IRAS-based Ltir is computed using a single poly- 
nomial equation (Eq[T]), the agreement shows it measured the IR 
LF very well. Perhaps, since a LF is an integrated quantity, it can 
be measured more reliably than Ltir of individual galaxies, even 
using data up to only lOO/xm. 

Our LF in Fig|51agrees well with lRush. Malkan and Spinogliol 

(Il993h . after correcting their results from a Hubble constant of 50 to 
75. Their LF slightly underestimate LF, perhaps because they only 
integrated IRAS photometry up to lOO/xm. Our LF also agrees well 
with Takeuchi, Yoshikawa, & Ishii ( 200 3), once their 60/im LF is 
converted to Ltir by multiplying 2.5 ( Takeuchi et al.i2 006). Their 
LF was measured in a SED model-free way, using different density 
estimators. It is reassuring that LFs measured in completely differ- 
ent ways agree well. See Sedgwick et al. (in prep.) for a similar 
attempt in the AKARI Deep field so uth. 

Following Is anders et al.l (l2003h . we fit an analytical function 
to the LFs. In the literature, IR LFs were fit better by a double- 
power law (Babbedge et al. 2006; Goto et al. 201 (1) or a double- 
exponential dSaunders et al.ll99G : IPozzi et al .120041 : [Takeuchi et al.l 



l2006l : iLe Floc'h et"aD l2005h than a Schechter function, which 
steeply declines at the high luminosity and underestimates the num- 
ber of bright galaxies. In this work, we fit the TIR LFs using a 
double-power law (Babbedge et al. 2006). 



^{L)dL/L* 



dL/L% {L<L* 



^{L)dL/L'' = 



1-/3 



dL/L\ [LyL"") 



(16) 



(17) 



o 10- 



10" 



• This work (CHEL model) 
X Sanders + 05 
+ Rush + 93 

9 10 11 12 

log [L,p(Lo)] 



Figure 6. Infrared luminosity function of the RBGS. The Ltir is mea- 
sured using the AKARI 9,18,65,90,140 and 160/im fluxes through an SED 
fit. Errors are computed using 1000 Monte Carlo simulations, added by 
Poisson error. The dott ed lines show the be st-fit double-power law. The 
crosses show data from lSanders et al.l (l2003h , who measured Lttr using 
IRAS photometry. The plus show data from iRush. Malkan and Spinogliol 
(Il993h . who measured Lpiji by integrating IRAS fluxes over 12-100/im. 




10 11 12 13 
log [LT,R(Le)] 



Figure 7. A luminosity histogram of galaxies used to compute Fig|6] 

Free parameters are L* (characteristic luminosity, L©), 0* (nor- 
malization, Mpc~^), a, and [3 (faint and bright end slopes), respec- 
tively. The best-fit values are summarized in Table |2] The local LF 
has a break atL*=7.8=b0.2 x 10^°. Understanding how this break 
(L*) evolves as function of cosmic time, and what causes the break 
is fundamental to galaxy evolution studies. This work provides an 
important benchmark in the local Universe. 

4.4 Bolometric IR luminosity density based on the TIR LF 

One of the primary purposes in computing IR LFs is to estimate the 
IR luminosity density, which in turn is a good estimator o f the dust- 
hidden cosmic star formation density (lKennicuttlll998h . provided 
the AGN contribution is removed. The bolometric IR luminosity of 
a galaxy is produced by thermal emission of its interstellar matter. 
In SF galaxies, the UV radiation produced by young stars heats 
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Table 2. Best double power-law fit parameters for the AKARI LFs, and the infrared luminosity density obtained from the fit 



Sample 




(/)*(Mpc-3dex-i) 


a (faint-end) 


/3 (bright-end) 




Total 


7.8±0.2 X IQio 


0.00037±0.00005 


1.8±0.1 


3.6±0.1 


8.5+^-J X 10^ 


SFG 
AGN 


10±0.2 X IQio 
7.8±0.2 X 10^0 


0.0001 8±0.00006 
0.00011 ±0.00005 


1.9±0.1 
1.7±0.1 


4.4ib0.3 
3.2±0.3 


8.0+2'5 X 10'^ 
l.6tr, X 10^ 



the interstellar dust, and the reprocessed light is emitted in the IR. 
For this reason, in star-forming galaxies (SFG), the bolometric IR 
luminosity is a good estimator of the current SFR (star formation 
rate) of the galaxy. 

Once we measured the LF, we can estimate the total infrared 
luminosity density by integrating the LF, weighted by the luminos- 
ity. We used the best-fit double-power law to integrate outside the 
luminosity range in which we have data, to obtain estimates of the 
total infrared luminosity density, Qtir- 

The resulting total luminosity density is Qtir= i^-^-l'l) x 
10^ LqMpc~^. Errors are estimated by varying the fit within 
1 cr of uncertainty in LFs. Out of Qtir, 7=bl% is produced by 
LIRG (Ltir > lO^^Lo), and only 0.4±0.1% is by ULIRG 
(Ltir > 10^^ Lq). a very small fraction of ^tir is produced by 
luminous infrared galaxies at z=0.0082, in stark contrast to high- 
redshift Universe. We found that ~ 30% of Qtir originates from 
Ltir < lO^^L©, where we do not have data and had to rely on 
the extrapolation of the faint-end tail of the LF. Therefore, simi- 
lar amount of uncertainty cannot be ruled out if the faint-end slope 
changes significantly under Ltir < lO^^L©. We will discuss the 
evolution of Qtir in Section [5^ 



5 INFRARED LUMINOSITY DENSITY AND ITS 
EVOLUTION 

5.1 Separating IR contributions from AGN and SFG 

We showed a local IR LF in Fig[6l However this includes both IR 
emission from star formation activity and from AGN. Since differ- 
ent physics govern star formation and AGN, it is essential to sep- 
arate their IR emission to understand the cosmic evolution of each 
component. For example, the cosmic star formation history cannot 
be addressed without subtracting the contribution from the AGN, 
and vise versa. 

Separating AGN from starbursts is not a trivial task. Many 
AGN/SF separation methods have been proposed such as X-ray, 
radio luminosity, optical line ratios, PAH strengths, submm proper- 
ties and so on, many of which often disagree with each other. We 
do not have a complete diagnosis of AGN/SB such as emission line 
ratios and X-ray for all of our sample. Therefore instead of classi- 
fying individual galaxies into AGN and SF, we attempt a statistical 
approach to separate contributions to LFs by AGN and star-forming 
galaxies. 

In this work, we simply us e a re sult of a new classifica- 
tion of I Yuan. Kewlev. & Sander si (l201(t) . wh ich used three care - 
fully examined optical line ratio diagrams (lKewlevetal.ll2006h 
to classify individual galaxies. Based on optical line ratios they 
classify IR galaxies into AGN (LINER, Seyfert 1 and 2), star- 
forming, and composite as a function of Ltir in their Ta- 
ble 4 and Figure 4. Their results agree well with other work 
(iTakeuchi. Yoshikawa. & !shiill2003l : [Gotoll2005h . within 20% or 
so. We use their fractions of AGN/SFG as a function of Ltir, 



to separate our LF into that of AGN and SFG. In this process, 
it is an open question how to handle composite galaxies, which 
are close to 50% of galaxies in almost all bins. In this work, we 
simply assigned the AGN/SF fractional ratio of that luminosity 
bin to separate composite galaxies (and ambiguous galaxies) into 
AGN/SF galaxies. We understand this is not an ideal approach. 
There will be up to ^50% uncertainty if these composite galax- 
ies are all AGN or SFG. It merely provides one method to sepa- 
rate IR LFs into those of AGN and SFG. A method to assess IR 
contribtuion by AGN and SFG is becoming realistic. For example, 
iGandhi et alJ too^ showed a tight correlation between nuclear 
Li2.3f^m and Lx-ray In the near future, more realistic classifi- 
cation of composite galaxies will become possible. We show the 
resulting AGN fraction from this process as a function of Ltir 
in FiglU At Ltir > lO^^L©, this procedure assigns 90% of IR 
luminosity to AGN. At Ltir < 10^° L©, the AGN fraction is 
20% or less. Th i s also agrees well with a recent measurement by 
iKartaltepe etaP feOlOl) . 

By applying this AGN/SF separation, in Fig|9] we show LFs 
separately for AGN and SFG. The total IR LF shown in Figl6l is 
also shown for comparison. As expected from Fig|8j Fig|9l shows 
that at Ltir > 10^^ L©, the AGN are responsible for almost all 
of the Ltir, forming a very steep bright-end drop for the SFG LF. 
On the other hand, at Ltir < 10^^ L©, SF galaxies explain most 
of the Ltir- 

Next, we fit a double-power law (Eqs.[T6]and[T7]) to the AGN 
and SF IR LFs, exactly as we did for the total IR LF in Figl6l The 
best-fit parameters are summarized in Table [3 The most notable 
difference is at the bright-end slope (/3), where the AGN LF has a 
very shallow slope of /3=3.2zb0.3, while the SFG has a steep slope 
of /3 =4.4±0.3. As expected, the faint-end tail (a) is not much 
different, in fact, being consistent with each other within la. The 
L* of SFG becomes brighter due to the steepening of the bright-end 
slope. 

The total infrared luminosities of AGN and SFG's were 
also m easured in separated LF's by iRush. Malkan and Spinogliol 
Adding their Seyfert 1 and Seyfert 2 LF's, and correcting 
for different Hubble constant, they get somewhat lower AGN/SFG 
LF ratios than our overall value of 1.6 : 8.0. Their methodology was 
simpler, and did not account for the population of LINER AGN in 
their 12/xm- selected galaxy sample. 

What do these differences in LFs bring to the IR luminosity 
density, ^ir, by AGN and SFG? We estimate the total infrared lu- 
minosity density by integrating the LFs weighted by the luminosity, 
separately for AGN and SFG. We used the double power law out- 
side the luminosity range in which we have data, to obtain estimates 
of the total infrared luminosity density, Q^tir, for AGN and SFG. 

The resulting total luminosity density {O^ir) is, 
llf|G=8.0l^:^ X 10" L0Mpc-^ and 1.6l°J x lO" 

LoMpc"^, as summarized in Table [3] Errors are estimated by 
varying the fit within la of uncertainty in LFs. These are also 
summarized in Table [3 The results show that among the total IR 
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Table 3. Local IR luminosity densities 



Sample 


n^^G (L0Mpc-3) 


nfg"^ (L0Mpc-3) 


Total 
LIRG 
ULIRG 


(8.0+2'5) X 10'^ 
(3.8l°-^) X 10^ 


(1.6^^3) X 10" 
(2.1+0;^) X 106 

(12+^) X 10^ 



luminosity density integrated over all the IR luminosity range, 
83% { ^AGN^oSFG ) of IR luminosity density is emitted by the 

IR IR 



, o^T^) is by at z=0.0082. The 

IR 



SFG, and only 16% ( 

results shows that at low redshift (z=0.0082), majority of the IR 
luminosity density is emitted by the SFG. Therefore, even if you 
convert all the Qir into SFR density, you will only overestimate by 
20%. However, the situation is different at higher redshift, where 
a number of papers reported luminosity evolution in the IR LF 
dPerez-Gonzalez et al.ll2005l : iLe Floc'h et"aDl2005l : iMagnelli et aP 
l2009h . The AGN/SFG separation will have much larger effect and 
thus more important at higher redshift. 

Once we have ^ir^, we can estimate star formation den- 
sity emitted in infrared light. The SFR and Ltir is related by the 
following equation for a Salpeter IMF, (m) oc m~^'^^ between 
0.1 - lOOM© (Kennicutt 1998). 



SFR[MQyr- 



1.72 X lO-^"" Ltir[Lq\ 



(18) 



By using this equation, we obtain SFR density = 1.3±0.2 

10-^MQyr-\ 

If we limit our integration to ULIRG luminosity range 
(Ltir > lO^^i^o), we obtain, Qff ^(ULIRG)= (1.5+i o) x 
10^ L0Mpc-^ and Q^^^(ULIRG)= (12^^) x 10^ LoMpc-^ 
In other words, at ULIRG luminosity range, AGN explain 



nf^'"' (ULIRG) 
nfGN (u LI RG)+nf^G (ULIRG) 



) of IR luminosity, again showing 



the AGN dominance at the bright-end. 
In the LIRG luminosity (Ltir 

X 10^ LoMpc 



Q^^^{URG)= (3.8 



+ 0.4X 
0.9; 



> 10 

3 



Lq), results are 
and Q^^^(LIRG)= 



(2.1+0 5) X 10^ LoMpc"^. This shows AGN contribution is al- 



-,AGN 
IIR. 



(LIRG) 



ready down to 35% ( nf^^ iULiRG)+nii^<^(LiRG) 
ity in the LIRG range. 



) of IR luminos- 



5.2 



Evolution of Qff ^ 



We have separated the ^jr^ from Qf^^ . Now we are ready to 
examine the evolution of Qsfr without contribution from AGN. 
In Fig[TOl we plot the evolution of as a function of redshift. 

Higher redshift results are taken from lGoto et al.l (1201 0), who also 
tried to exclude AGN using SED fitting to individual galaxies. Re- 
sults from the Spitzer survey and GALEX survey are also plotted. 
The ^jR^ shows a strong evolution as a function of redshift. The 
best-fit linear relation is Qir*^ (x(l-hz)^ °+°-^. This is consistent 
with most of earlier work. For example. |Le Floc'h et al.l (l2005l) ob- 
tained 7 = 3.9=b0.4 up to z^l.lPerez-Gonz alez et al.l(l2005h found 
7 = 4.0 ± .2 from z^O to 0.8. Babbedse ltaD (l2006h obtained 
7 = 4 5lnL MagnelU et al. (2009) obtained 7 = 3.6 =b 0.4 up to 
z=1.3.|^dighiero et al. (2010) found 7^8.8 ± 0.4 in the redshift 
interval of < z < l. lGruppioni et al.l(EoiO ) found -y 3.8 ±0.3 
up to z^l, with some evidence of flattening up to z^2. 

Once the IR luminosity density is separated into ULIRG and 
LIRG contribution, we found Q^^^ (ULIRG) (x(l-hz)^•^^°•^ 




10 11 
log [LT,R(Le)] 



12 



Figure 8. Fractions of AGN is shown as a function of Ltir. The ori ginal 
classification was taken from Fig.4 of Yua n. Kewlev. & SanderslbOlOh . and 
interpolated by a spline curve to be applied to our LFs. 
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Figure 9. IR LF is separated for star-forming galaxies (blue triangles) and 
AGN (red squares) using Fig|8] Total IR LF is shown with the black plus 
sign. The dotted-lines are the best-fit double-power laws. 



and n^^^{LIRG) oc(\+zf-^^^-^. Q^^^ (ULIRG) shows more 
rapid evolution than Qj^^ (LIRG), showing importance of lumi- 
nous IR sources at high redshift. 



5.3 Evolution of Qfg^ 

In turn, we can also investigate the evolution of 0.fR^ . This has 
been difficult in the literature since it is difficult to separate faint 
AGN individually in the presence of a host galaxy. By integrating 
IR LF AGN in FiglH we obtained Qfg^ =1.6^0.1 x 10^ LoMpc-^ 
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Figure 10. Evolution of TIR luminosity density by star-forming galaxies. Results from this work is plotted with triangles at z=0.0082. The red, blue and 
orange triangles show IR luminosity density from all galaxies, from LIRG only, and from ULIRG only. Higher redshift results in the solid lines are from the 
AKARI NEP deep field ( Goto et al. 2010). Shown with different colors are TIR luminosity density based on TIR LFs (red circles), S/xm LFs (stars), and llfim 
LPs (filled triangles). The blue open s quares and orange filled squares are for only L I RG and ULIRGs, also based o n our Lttr LFs. O verpl otted dot-dashed 
lines are estimates from the literature : lLe Floc'h et atl <2005h . lMagnem et al.l <2009h . |Perez-Gonzalez et al.l l2005h.lCaputi et al.l <2007h. an d lBabbedge et al] 
l2006h are in cyan, yellow, green, navy, and pink, resp ectively. The purple dash- dotted line shows UV estimate bv lSchiminovicl^t al.U2005h . The pink dashed 
line shows the total estimate of IR (TIR LF) and UV iSchiminovich et al.ll2005h . 



iKartaltepe et alJ (l201Qh recently reported that the mean SED 
of 70/im- selected sample from the COSMOS survey at z ~1 is 
similar to what has been observed locally. As a natural conse- 
quence, the fraction of AGNs as a function of Lir is in excellent 
agreement with that in the local universe (their Fig. 30). If so, we 
can apply the fAGN-LiR relation described in Fig|8]to IR LFs at 
higher redshifts to investigate the evolution of Qf^^ . 

Practically, we applied the fAG N-Lin relation in Fig|8]to IR 
LFs at < z < L5 p resented bv iGoto et al. (2010) to obtain 
Q^g^. lGotoetaPbOlOh used their own AGN/star-forming galaxy 
classification to remove AGN from their LFs, but for a fair com- 
parion, we applied the same methodology as used in Section 15.11 
to their LFs before they subtract the contribuion from individual 
AGN. 

In Fig[TT] we show the evolution of Qf§^ , which shows a 
strong evolution with increasing redshift. At a first glance, both 



^fn^ and ^jr^ show rapid evolution, suggesting that the cor- 
relation between star formation and black hole accretion rate con- 
tinues to hold at higher redshifts, i.e., galaxies and black holes 
seem to be evolving hand in hand. However, there are some dif- 
ferences as well. When we fit the evolution with (1-i-z)^, we found 
(x(l-i-z)^"^^°"^, which is possibly more rapid evolution than 
^iR^ although errors are significant. The contribution by ULIRGs 
quickly increases toward higher redshift; By z=1.5, it exceeds that 
from LIRGs. Indeed, we found Qfg^ (ULIRG) (x(l-hz)^^^°^ 
and Qf^^'iLIRG) (x(l-hzf -^^^ ^ 

It seems there is no sign of flattening at high redshift 
in Qf^^{z). It would be interesting to c ompare this with a 
number-density evo lution of optical QSOs dCroom et al.l l2004l : 
[Richards et al.ll2006h. which peaks at z=2-3, and the evolution of 



X-ray AGN dUeda et aLbOOllHasinger. Mivaii. & Schmidtl2005h . 
lAird etalJ tOl& i recently presented the evolution of X-ray lu- 
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Figure 12. X-ray luminosity density to ^f^^ ratio is plotted as a func- 
tion of redshift. The ratio is normalized at the minimum value. The X-ray 
luminosity densities of black circles and blue square s are f rom lAird et al.l 
(l201Qh (2-10 keV) and iHasinger. Mivaii. & Schmidt! (l2005h (0.5-2 keV) , 
respectively. 



minosity density, which shows a turnover at around z~1.2zb0.1. 
In FigdU our ^fR^ does not show any sign o f decline at 
least u p to z=1.5, in contrast to the X-ray results. iBarger et"aD 
(l2005h also showed much shallower evolution in X-ray than our 
IR results. In FigIT2l we investigate the evolution of X-ray lu- 
minosity de nsity to r a tio. X-ray l umino sity densities are 



taken from lAird et all dioTol) lAird et all (l201Qh (2-10 keV) and 
iHas inger. Mi vaji. & Schmidt! tOOS h (0.5-2 keV). It is interesting 
that the ratio is consistent with a constant value at < z < 1, and 
shows a possible increase at 2; > 1, although it is unfortunate our 
errors are too large to draw a firm conclusion. The possible differ- 
ence may suggest an increase of obscured AGN light toward higher 
redshift at 2; > 1 compared with optical/X-ray unobscured AGN. 
Post-starburst AGN galaxies found by Goto (2006) may be an ex- 
ample of such phenomena. Recently, iTreister et al.l (I2OIOI) showed 
an increase in the number ratio of obscured/unobscured QSOs with 
increasing redshift, concluding that a merger-driven black hole evo- 
lution model is consistent with the observed result. Our Fig|T2lmav 
suggest a similar obscured/unobscured AGN evolution in terms of 
AGN luminosity density, instead of the number density. A possible 
deviation from ^jr^ also suggests AGN may have formed earlier 
than star-formation in terms of infrared luminosity density. More 
accurate comparison of X-ray and IR luminosity density evolution 
of AGN will bring an interesting physical implication on the sub- 
ject. In particular, reducing measurement errors at high redshift is 
urgent. 



6 SUMMARY 

Using AKARI's 6-band IR photometry in 9, 18, 65, 90, 140, and 
160/im with much improved spatial resolution to resolve back- 
ground cirrus/confusing sources, we have re-measured Ltir of the 
RBGS through IR SED model fitting. This is a significant progress 



considering that the IRAS used one linear equation to estimate 
Ltir, with only data at < lOO/xm. 

By using this new Ltir measurement, we constructed local 
IR LFs separately for SFG and AGN. We also computed local in- 
frared luminosity density through the derived LFs, and compared 
^iR^ and n to those at higher redshifts. 

Our findings are as follows. 

• SED model-to-model variation in estimating Ltir is less than 
25%. 

• We present Lmir-^o-Ltir conversions for Spitzer 8/im, 
AKARI 9/im, IRAS 12/im, WISE 12/im, ISO 15/im, AKARI 
18/im, WISE 22/im, and Spitzer 24/im filters. These conversions 
provides us with a useful tool to estimate Ltir with a MIR band 
only. 

• Re-constructed local TIR LF with the AKARI data is consis- 
tent with that from the IRAS, i.e., the AKARI's better data show 
the IRAS measurement was correct. 

• By integrating the IR LF weighted by Ltir, we obtain the 
local cosmic IR luminosity density of VLtir= {^.5t}2%) x 10^ 
LoMpc"^. 

• LIRG and ULIRG contribute to VLtir a little; Only l±l% 
of VLtir is produced by LIRG {Ltir > 10^^ Lq), and only 
0.4±0.1% is by ULIRG (Ltir > lO^^L©) in the local Universe, 
in stark contrast to high redshift results. 

• Compared with high redshift results from the AKARI 
NEP deep survey, we observed a strong evolution of 



^rm oc(l-hz) 



4.0±0.5 



, after removing AGN contribution. 



• We showed an evolution of ^fR^ scales as oc(l-i-z)^'^^°'^. 
ULIRG contribution to the Hj^^ exceeds that by LIRG by z=1.5. 

• ^fR^/^jfi^ay ratio shows a possible increase at z > 1. If 
confirmed, this may suggest an increase of obscured AGN at 2; > 1. 
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Figure 11. Evolution of TIR luminosity density by AGN. Results from this work is plotted with stars at z=0.0082. The red, blue and orange poin ts show IR 
luminosity density from all AGN, from LIRG AGN only, and from ULIRG AGN only. Higher redshift results are from th e AKARI NEP de ep field JGoto et alJ 
bOlOl) . with contribution from star forming galaxies removed. Brown triangles are ^f^^ computed from the 8/im LFs ( iGoto et al.ll201Ql) . 
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